Enhancement of the v e flux from astrophysical sources by two photon annihilation 

interactions 
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The ratio of anti-electron to total neutrino flux, <& Pll : $„, expected from jry interactions in 
astrophysical sources is < 1 : 15. We point out that this ratio is enhanced by the decay of 
pairs, created by the annihilation of secondary high energy photons from the decay of the neutral 
pions produced in py interactions. We show that, under certain conditions, the <J>£; E : <E>„ ratio may 
be significantly enhanced in gamma-ray burst (GRB) fireballs, and that detection at the Glashow 
resonance of i> e in kilometer scale neutrino detectors may constrain GRB fireball model parameters, 
such as the magnetic field and energy dissipation radius. 

PACS numbers: 96.40.Tv, 14.60.Pq, 98.70.Rz, 98.70.Sa 



I. INTRODUCTION 

High energy neutrinos are expected to be produced in 
astrophysical sources mainly by pry interactions, leading 
to the production and subsequent decay of charged pi- 
ons: 7r + — > e^vjj^v^ (see, e.g., 1] for recent reviews). 
Neutrino oscillations lead in this case to an observed ra- 
tio of v e flux to the total v flux of ~ 1 : 15 Q (or lower, 
in case muons suffer significant electromagnetic energy 
loss prior to decay |3]). For neutrinos produced in in- 
elastic pp (pn) nuclear collisions, where both 7r + 's and 
7r~'s are produced, the ratio is ~ 1 : 6, and it was sug- 
gested that measurements of the v e to v e flux ratio at 
the M^-resonance may allow one to probe the physics of 
the sources by discriminating between the two primary 
modes of pion production, pj and pp collisions A] . This 
test for discriminating between the two mechanisms is 
complicated by the fact that the ratio of $ Pe to <£>„ pro- 
duced in pry interactions can be enhanced to a value sim- 
ilar to that due to inelastic nuclear collisions in sources 
where the optical depth to py interactions is large (e.g. 
0). In this case, neutrons produced in pry — > nn + inter- 
actions are likely to interact with photons and produce 
tt~ before escaping the source, leading to production of 
roughly equal numbers of 7r + 's and 7r~'s. 

In this paper, we point out that the : <!>„ ratio 
from pry interactions may be enhanced above 1:15 also 
in sources with small py optical depth. Neutral pions, 
which are created at roughly the same rate as charged 
pions in pj interactions, decay to produce high energy 
7-rays. These 7-rays typically carry ~ 10% of the initial 
proton energy, and may therefore interact with the low 
energy photons (with which the protons interact to pro- 
duce pions) to produce /i + pT pairs. The decay of muons 
yields (after vacuum oscillations) $p c : <!>„ ~ 1 : 5, thus 
enhancing the v e fraction. 

We discuss below a specific example, the widely consid- 
ered fireball model of GRBs. In this model, the observed 
7-rays are produced by synchrotron radiation of shock 
accelerated electrons in the magnetic field which is as- 



sumed to be a fraction of the total energy (see e.g. 
for reviews). The protons are expected to co-accelerate 
with electrons to ultra- high energy , and produce high 
energy neutrinos by pj interactions 0. We calculate 
below the additional neutrino flux, due to the decay of 
muons produced by secondary photon annihilation, for 
a typical long duration GRB, and show that the en- 
hanced D e flux may be detectable at the Glashow res- 
onance (P e e — > W~ — > anything 8J) in kilometer scale 
neutrino detectors such as IceCube |!j 

The enhancement of <£> Pe : <!>,, due to 77 interac- 
tions in pry sources makes the discrimination between pj 
and pp neutrino sources more difficult. On the other 
hand, it may provide a new handle on the physics of the 
source. We show below that for GRBs the enhancement 
of D e flux depends on model parameters which are poorly 
constrained by observations, namely the magnetic field 
strength and the energy dissipation radius. Detection of 
u e 's at the Glashow resonance, in conjunction with 7-ray 
detection, may therefore constrain these parameters. 



II. FIREBALL MODEL AND PHOTON 
SPECTRUM 

The minimum observed GRB fireball radius r may 
be estimated by requiring that it is optically thin to 
Thomson scatterings: r^ h = a^n'r' < 1 (denoting 
the comoving and local lab. frame variables with and 
without a prime respectively). Here n' is the density 
of scatterers in the fireball, r' — r/T is the size of 
the interaction region and Y is the bulk Lorentz fac- 
tor. The radius at which r^ h w 1 is the photospheric 
radius r p h. For a kinetic luminosity of the fireball, 
mostly carried by the baryons, the number density of 
the baryons, and of the leptons which are coupled to 
the baryons, is n' b « Lk/(4-7rr 2 r 2 m p c 3 ). The observed 
isotropic equivalent 7-ray luminosity of a long duration 
GRB is L 52 = L 7 /10 52 erg/s ~ 1. Assuming L 1 = e e L^ 
with e e r-j 0.05e ej _i.3 (a parametrization which is moti- 
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vated below), the photospheric radius is 
a Th L 7 /e e _ 7 4 x 1Ql2 L 52 



47rr 3 m„c 3 



^,-1.3^.5 



(1) 



for T2.5 = r/316 ~ 1. The radius at which the bulk 
kinetic energy dissipation occurs, e.g. by internal shocks, 
is in general r > r p h- 

The 7-ray spectrum of a GRB fireball at a dissipation 
radius r = 10 14 ri4 cm peaks at a typical energy 

e 7 , P k = HcT 2 (3j 2 min qB')/(2m e c 2 ) 

~ 500 (e^ 1 . 3es ^ lJ L 52 ri 5 /r2 4 ) 1/2 keV, (2) 

due to synchrotron radiation by electrons with a Lorentz 
factor 7e min ~ e e (m p /m e ). Here, "f' emin is at the lower 

end of a a l/7 e p distribution of electron Lorentz fac- 
tor, with p > 2, created by Fermi acceleration in the 
shock. The magnetic field is assumed to be B 2 /8tt « 
e_B-f'k/(47rr 2 r 2 c), where Eb ~ 0.1£_b^_i is the equipartion 
value, currently unconstrained in the GRB prompt phase. 
Note that e 7 , p k oc 1/r will be larger than the above value 
for r w r p h, with other parameters fixed. 

For a GRB at a luminosity distance cLl the observed 
7-ray spectrum is generally approximated with a broken 
power-law Band fit ^} , 



de 7 47rd|e 2 , pk 1 (e 7 /e 7 , pk ) 2 ; e 7 > e 7 , P k- 



(3) 



The spectrum deviates from this at low energy, becoming 
dN 1 /d e 7 oc e 7 ^ for e 7 < e 7 , S ai the energy below which 
synchrotron self-absorption becomes dominant. Theoret- 
ical modeling indicates a value 
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2-4 r 2 7 ' in <r[gfic] 4 B 2 /[m 3 e c e ] 



1/3 



( £ B,-iii 2 /[ee,-i. 3 rl 5 r? 4 ]) 1/3 keV ) (4) 



for p = 2. The differential number density of photons is 



dN^/de'^ « L 7 /(47rr 2 ce 2 !pk ) 



( £ 7 ,sa/ e 7:P k) ( e 7 / £ ' 
( £ 7 / £ 7 ,pk) 
^ ( £ 7 / £ 7 ,pk) 



,3/2 



7,sa y 



■ e' < e' 



7,pk ^ 7 ^ 7,sa 
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Electron synchrotron radiation produces a power law 
7-ray spectrum at energies above e 7 , p k [see Eq. JjSj] 
which depends on the maximum Lorentz factor. Other 
mechanisms can contribute to an extension of the 7-ray 
spectrum in Eq. iJSJ to high energies. High energy elec- 
trons can inverse Compton scatter synchrotron photons 
up to an energy similar to the maximum shock accel- 
erated electron energy (which we derive shortly) in the 
Klein-Nishina limit in one mechanism. Here we consider 



ultra-high energy 7-rays from w° decays which are pro- 
duced by 7)7 interactions of shock accelerated protons 
with synchrotron photons as pj — » A + — * pir° — > P77. 

The maximum proton energy is calculated by equating 
its acceleration time i^cc ~ e p/(9 c -S') to the shorter of 
the dynamic time t' dyn w r/(2Tc) and the synchrotron 

cooling time t' » QirrripC 3 / {a^m^e^B 2 ) as 



e P ,max = (&nm p c qT 2 /[<j Th m 2 e B']) 

« 3.3 x 10 11 ( £e ,_ 1 . 3 r6 5 r 2 4 /[e B ,_ 1 L 52 ]) 1/4 GeV 



qB'r 



5.5 x 10 1 



11 ( Eg, -1L52 
£ e,-ir 2 5 



1/2 



GeV, (6) 



respectively for t' acc = t' syn and t' acc = t' dyn . For elec- 
trons, t' dyn 3> t' syn typically and the maximum elec- 
tron energy, using Eq. © for electrons, is e e ,max ~ 



ltJ E e -1.3^,-1^52 i 2.5 r 14 UeV - 

At a given incident proton energy e p , the threshold 
photon energy leading to a A + resonance interaction is 



-1/4 r -l/4 F 3/2 1/2 



C 7 A + 



o.3r 2 



(ep/GeV) 



GeV. 



(7) 



The optical depth for this interaction may be calculated 
using a delta- function approximation with a cross-section 
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7 ,A+ 
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using Eq. (J5J. Note that the target photon spectrum, 
within parenthesis, is now evaluated at the A + resonance 
energy [see Eq. Q)] for an incident proton energy e' p . In 
particular, Eq. J7J) may be used to replace A+ by e' p 
in Eq. (|5j). As a result, the optical depth in Eq. (JSJ is 
expressed as a function of e' . The spectral shape of the 
optical depth is then oc (e p ) 9_1 , where q is the spectral 
index, dN^/de'^ oc (e 7 )~ 9 , in Eq. ©. Also the order 
is reversed, i.e. q = 2, 1, —3/2 in r p7 oc (e p ) 9_1 instead 
of q = —3/2, 1, 2 in dN' /de^ oc (e 7 )~ 9 , according to the 
condition in Eq. (7J|. 

Protons lose ~ 20% of their energy by P7 interactions 
to 7T° and e 7 w 0.1e p for each secondary photon. With 
an equal probability to produce 7r° and 7r + in each pj 
interaction, the resulting ir° decay photon flux is 



— = min[iy P7 ; 



0.2 (^/e e )i 7 



(9) 



Here £ p is the proton fraction undergoing shock acceler- 
ation. For £ p = 1 and e e = 0.05 we have r^ 7 rs 1 which 
leads to the observed flux level in Eq. Secondary 
pions from A + decay, and subsequent decay photons and 
neutrinos follow the dN p /de p oc e~ p spectral shape of the 
protons for a constant optical depth. For an optical depth 
of spectral shape oc e^" 1 , the resulting pion, photon and 
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neutrino spectra would be dN/de a e g 1 p . The w° de- 
cay photon spectrum in Eq. © is then dN 1 /de 1 oc e~ 2 
between e 7 = 0.03r 2 /e 7 , pk GeV 2 and 0.03r 2 /e 7 , sa GeV 2 , 
oc e" 1 below e 7 = 0.03r 2 /e 7iP k GeV 2 and oc e 7 9 ^ 2 above 
e 7 = 0.03r 2 /e 7iSa GeV 2 due to self-absorption following 
Eqs. 10 and ©. 

Note that, the luminosity of shock- accelerated protons 
is l/e e = 20 times the shock-accelerated electron lumi- 
nosity. In the fast cooling scenario, valid in the GRB in- 
ternal shocks, the electrons synchrotron radiate all their 
energy into observed 7-rays. Thus L p L 1 /e e . In 
a single interaction the secondary pion (charged or 
neutral) luminosity is L w w 0.2L p « 4L 7 /e ei _i.3. The 
neutrino luminosity of all flavors from charged pion de- 
cay, assuming equal energy for all 4 final leptons, is 
L v « (l/2)(3/4)L 7r w 1.5L 7 /£ e _i. 3 . The 1/2 factor 
arises from the equal probability of tt° and tt + produc- 
tion. The neutrinos carry away energy from the fireball, 
and the rest of the pion decay (e + from tt + and 77 from 
7r ) energy is electromagnetic (e.m.), with a luminosity 
L c . m . w L n — L v w 2.5i 7 /e e ,-i.3- A significant fraction 
of the 7r°-decay L 7 [Eq. 0] would be converted to muon 
pairs and subsequently to neutrinos, as we discuss next. 
A substantial (small) fraction of the rest of £ c .m. would 
be emitted at r > r p h (r ~ r p h) as low energy photons 
with a luminosity not significantly above the observed 
7-ray luminosity. These, however, do not affect substan- 
tially the neutrino flux calculated from very high energy 
7-rays interacting with soft photons. 



tion opacities by integrating Eq. (|10fl piecewise as 
T 77 ^ ;+i -(e 7 ) = r 2 m 2 c 3 L 7 /(8re 2 pk e 7 ) 

1 t <W, P k A 2 _ ( e 7 , P k V 

2 \ ei.th / V £ t / 

C/,th 



m^c 4 r 2 
' £ 7 < 2e T , pk 



X < 



- 1 



> £ 7 > 2c 
3/2 

2 



m, c r 



(11) 
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;e 7 > 



2c- 



Here we defined the threshold energy for lepton pair pro- 
duction as e^th = m 2 c 4 r 2 /2e 7 . Note that the high energy 
photons produce e + e~ pairs dominantly at lower energy. 
The ratio of the two opacities k = t^^u+u,- / T 77 — >-e+e- 
becomes unity for higher energy photons, since the cross- 
section is the same for and e + e~ pair productions 
above the muon pair production threshold energy. 

A. Muon decay neutrino flux 

Muon pairs decay to neutrinos as /i~ — > e~v e v p and 
/i + — » e^Vepp shortly after they are created in the 
cm. frame of the 77 collision. In the observer's frame 

w e 7 /2, and the particle pairs move radially along the 
incident photon's direction. For simplicity we assume 
that the v e and created from //-decay carry 1/3 of 
the muon energy each. The observed neutrino energies 
are then e 7 /6 for each flavor. The neutrino source flux, 
which is the same for v ei , v e and previous to any 
flavor oscillation in vacuum, is 



III. TWO PHOTON PAIR PRODUCTION 



^ s __ dN v 

€v ^ " eu de v 



= miti[l,T 77 _ >/1 + M -]/c& 



dN 1 
de-y 



(12) 



High energy 7-rays can produce lepton pairs, 
{I = e,/i), with other photons which are above a thresh- 
old energy wth = mic 2 in the center of mass (cm.) 
frame of interaction. For an incident (target) photon 
of energy e' 7 ; (e 7 t ) in the comoving GRB fireball frame, 
lu = (2e 7ji e^ t ) 1 / 2 , and the cross-section for l + l~ pair pro- 
duction may be written, ignoring the logarithmic rise fac- 
tor at high energy, as tr 77 ^/+/- « irr 2 (mic 2 /uj) 2 , where 
r e is the classical electron radius. The corresponding op- 
tical depth is 



T 77( £ 7 a) = "F / a rv(^ e 7,t)-j^ de ' 7 ,t- (1°) 



7.t 



Given the power-law dependence of the photon distribu- 
tion in Eq. (JSJ, we may calculate the pair produc- 



The high energy muons produced from 77 interactions 
may lose a significant fraction of their energy by syn- 
chrotron radiation before they decay into neutrinos (with 
a decay time tdoc), if their energy is above a break energy 



^,sb 



= (6TTmlc 5 r 2 /[t dcc a Th m 2 e B'} 



,1/2 



5 x 10 7 (ee^Ti.rlJieB^L^ 2 GeV.(13) 



The corresponding neutrino break energy from muon de- 
cay is e„ iS b = ep.sb/6. For e„ > e v s b, the neutrino flux 
index would steepen by a factor 2 [lTj ■ 

We have plotted in Fig. the v e flux at the source, 
e 2 ^1 (same for v e , v e and P M ), previous to any 
vacuum oscillation, arising from 77 — > interac- 
tions and the associated muon decays, for a GRB of 
isotropic equivalent luminosity L 7 = 10 52 erg/s, which 
is average for a long GRB @. We assume a redshift of 
z ~ 0.1, which is near the low end of observed redshifts; 
there have been a few spectroscopic redshifts observed 
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in the 0.1 — 0.2 in the past 8 years, and indirect red- 
shift measures, such as time lags, indicate many more 
bursts in this range among the BATSE sample (see, e.g., 
Ref. U and references therein). Also plotted are the 
v e source flux: <? v <&% = min[l, Tp 7 ](0.2/8)i 7 /(47rd| i e e ), 
from j»7 — > A + — > mt + interactions and subsequent 7r + 
and /x + decays. Different panels are for different bulk 
Lorentz factor T and dissipation radii r. 




6 7 8 56789 10 

Neutrino energy log £ v (GeV) 

FIG. 1: Source flux of z/ e (same for u e and P M ) from 
77 — > interaction and subsequent /i-decays, compared 

to the canonical u e flux (same for v^) from single p7 — > nir + 
interactions and subsequent -k + , fj, + decays. This is for a long 
duration GRB of L 7 = 10 52 erg/s at redshift z ~ 0.1. The 
solid, dashed, dot-dashed and dotted lines are for magnetic 



field parameters eb = 10 ,10 , 10 1 
tively, for different F and r combinations. 



and 10 respec- 



B. Neutrino flavor oscillation and flux on Earth 

While neutrinos are created via weak interactions as 
flavor eigenstates, their propagation is determined by the 
mass eigenstates. The flavor eigenstates v a and the mass 
eigenstates z/j are mixed through a unitary matrix defined 
as v a = J2j U*jVj, where a = e, /i, r and j = 1,2,3 for 
three known flavors. The probability for flavor change 
by vacuum oscillation is given by V Va ^ Vfi — J2j Wi3j\ • 
U a j I , for the neutrino propagation from their sources 
to Earth over astrophysical distances. 

We use the standard expression for U a ,j with solar mix- 
ing angle 9q = 9\i = 32.5° and atmospheric mixing angle 
(9 at m = 023 — 45° |13|. The unknown mixing angle #13 
and the CP violating phase may be assumed to be zero 
given the current upperbounds from reactor experiments. 
Using these values for U a j and V Va ^ v „ results in a rela- 
tionship between the source neutrino fluxes <5> s u and the 
expected neutrino fluxes on Earth <f>„ which is given by 
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(14) 









For antineutrinos 7 3 p Qf _,p /3 is the same as above. 

Different production mechanisms produce v and v 
fluxes at the source with different flavor proportions. 
Their production ratios may be expressed as normalized 
vectors, shown in the left hand side of Eqs. I|15l fc Hul l. 
The corresponding flux ratios at Earth, using Eq. i|14fl . 
are shown in the right hand side of Eqs. (|15l fc lTrj|) below 
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The source ^-fluxes plotted in Fig. ^wiU be modified ac- 
cordingly Note that &v c {ll)/^vSPl) ~ 4 for the same 
initial 77 and fry flux levels. The i/ e -flux component is 
1/5 (1/15) of the total ^-flux from 77 (pj). The observed 
ratio of v e to total v fluxes from both the single 737 inter- 
actions and 77 interactions can be calculated using Eqs. 
(fH)l & Hrj|) from the source fluxes as 



£^ = 0.2$| e;p7 + 0.8$| e:77 



(17) 



We have plotted this ratio in Fig. [21 for different GRB 
model parameters. The ratio is enhanced from the canon- 
ical (pry) 1/15 value in the energy range where 77 interac- 
tions contribute significantly (see Fig. QJ. Interestingly, 
the enhancement takes place over a small energy range 
which may be explored to learn about the GRB model 
parameters as we discuss next. 



IV. NEUTRINO DETECTION 

We consider here the anti-electron neutrino detec- 
tion channel at the Glashow resonance energy e^rcs = 



m 



V2m e 

.3 



6.4 PeV 



The number of electrons in 
the 2 km a effective IceCube volume is iV e!C ff w 6 x 10 38 
and the corresponding number of downgoing v e events 
from a point source of flux is 



AtNp 



Am P c 2 



(18) 



Here g 2 ~ 0.43 from the standard model of electro-weak 
theory, and At is the duration of the emission. 

We have plotted in Fig. [21 the expected number of v e 
events at e^ros from a GRB fireball for various value of r, 
Eb and T. We have assumed here that shock accelerated 
protons interact once with synchrotron photons (two top 
panels), losing ~ 20% of their energy. For very high pr/ 
opacity the protons can lose most of their energy through 
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FIG. 2: Observed ratio of electron antineutrino flux to the 
total neutrino flux <£„ C /<1?„ for the py and 77 source fluxes 
plotted in Fig. This is for a long duration GRB of 
Lj = 10 52 erg/s at redshift z ~ 0.1, the solid, dashed, dot- 
dashed and dotted lines indicate magnetic field parameters 
eb = 10 _1 , 10 -2 , 10 -3 and 10 -4 respectively, for different 
r and r combinations. Note that the 1/15 ratio from single 
py — > nn + interactions is enhanced by 77 — > interac- 
tions in certain energy ranges which depend upon the GRB 
model parameters. 
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FIG. 3: Resonant v e events in IceCube from 77 and py inter- 
actions in a GRB, as a function of the magnetization sb for 
Lorentz factor F of 100 (solid lines), 178 (dashed lines), 316 
(dot-dashed lines), 500 (dotted lines) and fireball radii equal 
to the photospheric radius (two left panels) and ten times the 
photospheric radii (two right panels). We used single py in- 
teractions (two top panels) and multiple py and subsequent 
ny interactions proportional to r vl (two bottom panels) for 
comparison. Other GRB parameters are z — 0.1, At = 30 s, 
L 7 = 10 52 erg/s and e e = 0.05. 



the py and ny interaction chains. This could lead to 
the i^-fluxes from both py and 77 plotted in Fig. ^ to 
roughly increase by a factor five. The py source ^-fluxes 
may reach ratios <J>^ = <J>! = [1,2,0], from ir + and tt~ 
decays, which at Earth would be = $p = [1,1,1] 
for e„ < e^ sb/2. The resonant P e -events in such case, 
from both 77 and py fluxes, are plotted in Fig. [21 (two 
bottom panels). For simplicity, we assumed <£v oc $ 7 oc 
min[l, 0.2Tp 7 ] for the py z^-flux and 7-flux in Eq. I)12|l. 
Also $ p = [1, 1, 1] from p/ny interactions for t' p1 > 2. 

The background for astrophysical v e detection is 
mostly due to atmospheric prompt neutrinos from 
cosmic-ray generated charm meson decays. A 
parametrization for the z/ e and T> e atmospheric flux is 
given by (lif 



1.5x10" 



(jjatm 



;e v < 1.2 x 10 6 GeV 



l+3xl0- 8 e„ 

x GeV -1 cm" 2 s" 1 sr^ 1 . (19) 



The corresponding background z/ e -events at the Glashow 
resonance energy is < 10~ 7 for a GRB within a ~ 100 
s time window, allowing the full directional uncertainty 
(2tt sr), given the poor current knowledge of the v e or 
v e signal reconstruction in neutrino Cherenkov detectors, 
fromEq. flg) . 



observed $ Pe : $„ flux ratio. A different source of en- 
hancement of the v e flux may be pp interactions 13, 0] . 
In GRBs, however, the optical depth to pp is low 0, ex- 
cept in buried jets leading to v precursors 01 > where pp 
interactions are expected to lead to v's at energies ~ TeV. 
The number of resonant v e events arising from py inter- 
actions is essentially independent of eb (for £b < 10~ 2 ) 
for any T. On the other hand, the number of resonant 
v e events arising from 77 interactions varies significantly 
with r and r. It may become as large as the py con- 
tribution for 10~ 2 < £ B < 10- 3 , 100 < T < 300 and 
r ph r ?5 3r p h- For long bursts of average isotropic- 
equivalent luminosity at a redshift ~ 0.1, which from 
past experience are electromagnetically detected every 
few years, IceCube could probe the D e enhancement, and 
thus the value of the magnetization parameter and dis- 
sipation radius, by measuring the $ Pe : $> u flux ratio. 
Finally, we note that a moderate excess of v e events com- 
pared to and v T events may also be an indication for 
the presence of a 77 component. 



DISCUSSION 
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